Draft version February 1, 2008 

Preprint typeset using LAT^X style emulatcapj v. 04/03/99 



THE MERGER HISTORY OF SUPERMASSIVE BLACK HOLES IN GALAXIES 

Kristen Menou, 1 Zoltan Haiman, 2 and Vijay K. Narayanan 

Princeton University, Department of Astrophysical Sciences, Princeton NJ 08544, USA, 
kristen, zoltan, vijay@astro.princeton.edu 
Draft version February 1, 2008 



o 
o 



o 

(N 
> 

On 



ABSTRACT 

The ubiquity of supermassive black holes (SMBHs) at the centers of nearby luminous galaxies can 
arise from the multiple mergers experienced by dark matter halos in hierarchical structure formation 
models, even if only a small fraction of these galaxies harbor SMBHs at high redshifts. We illustrate 
this possibility using cosmological Monte Carlo simulations of the merger history of dark matter halos 
and their associated SMBHs. In our most extreme models, in order to populate nearly every bright 
galaxy with a SMBH at z = 0, only a few percent of the halos with virial temperatures above 10 4 K 
are required to harbor a SMBH at high redshift. This possibility must be included in studies of the 
luminosity function and the clustering properties of quasars. We predict the number of SMBH merger 
events that are detectable by the gravitational wave experiment LISA, as a function of redshift, out to 
z — 5. Although the event rates can be significantly reduced in scenarios with rare SMBHs, a minimum 
of ~ 10 detectable merger events per year is predicted if SMBH binaries coalesce efficiently. The observed 
distribution of events with redshift could yield valuable information on the SMBH formation process. If 
SMBH binaries do not coalesce, we find that at least several SMBH slingshot ejections probably occurred 
from z ~ 5 to the present in each galaxy more massive than ~ 1O 11 M at z — 0. Although our results are 
sensitive to the minimum cooling mass assumed for the formation of SMBHs, we expect the qualitative 
predictions of our models to be robust. 

Subject headings: black hole physics - cosmology: theory - quasars: general - galaxies: active, nuclei, 
interactions - gravitation - relativity 
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1. INTRODUCTION 

There is growing observational evidence for the pres- 
ence of supermassive black holes (SMBHs) at the centers 
of many, if not all, nearby, weakly-active or inactive galax- 
ies (Kormendy & Richstone 1995 and references therein; 
Magorrian et al. 1998; van der Marel 1999). Early at- 
tempts to relate the properties of the nuclear SMBH to 
those of the parent galaxy showed the existence of a lin- 
ear relation between the SMBH mass and the galaxy lu- 
minosity (or that of its bulge component in the case of 
disk galaxies; Kormendy & Richstone 1995; Magorrian et 
al. 1998), albeit with a large scatter (of about 0.5 dex in 
the ratio log[MsMBii/-Mbuige])- Recently, a substantially 
tighter correlation between the SMBH mass and the stellar 
velocity dispersion of the host galaxy bulge was reported 
by Ferrarese & Merritt (2000) and Gebhardt et al. (2000), 
suggesting that the formation and evolution of SMBHs and 
the bulge of the parent galaxy may be closely related. 

The nuclear SMBHs are most probably "dead quasars" , 
i.e., relics of quasar activity that may have occurred in 
many galaxies over their past history (e.g. Lyndcn-Bcll 
1967; Soltan 1982; Rees 1990; Richstone et al. 1998). At 
cosmological distances, the presence of SMBHs in galaxies 
is inferred from the powerful quasar activity phenomenon, 
which is commonly attributed to the radiatively efficient 
accretion of gas onto a SMBH of mass ~ 10 7 — 10 9 M 
(Salpeter 1964; Zel'dovich 1964; Lynden-Bell 1969; Rees 
1984; Blandford 1992). The simplest assumption relating 
the observations of nuclear black holes in both nearby and 
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high redshift objects is that a SMBH is present in every 
galaxy at any redshift, while only a fraction of these galax- 
ies are active at any given time. A number of models based 
on this assumption have been constructed over the years 
to explain the cosmological evolution of the quasar popu- 
lation with redshift (e.g. Efstathiou & Rees 1988; Small & 
Blandford 1992; Haehnelt & Rees 1993; Haehnelt, Natara- 
jan & Rees 1998; Haiman & Loeb 1998; Haiman & Menou 
2000; Kauffmann & Haehnelt 2000; Cavaliere & Vittorini 
2000) and were shown to be consistent with current obser- 
vational constraints at both low and high redshift. 

In this paper, we show that models in which only a small 
fraction of galaxies at high redshift harbor SMBHs are 
also consistent with the presence of SMBHs in the bulges 
of nearly all local, massive galaxies. This possibility arises 
because in standard hierarchical structure formation mod- 
els, massive halos experience multiple mergers during their 
lifetimes. This implies that the mass of central SMBHs at 
redshift z — could originate from just a few selected pro- 
genitor halos, rather than having been distributed more 
uniformally in all the progenitors. Low-frequency gravita- 
tional waves, emitted when the SMBHs at the centers of 
merged halos coalesce, are detectable to high redshifts by 
the planned Laser Interferometer Space Antenna (LISA; 
Folkner 1998). We estimate the event rates detectable 
by LISA in scenarios with rare galactic SMBHs at high 
redshift and we suggest that LISA could prove a useful 
tool in constraining the characteristics of the population 
of SMBHs over a range of redshifts. 

The outline of this paper is as follows. In §2, we briefly 
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describe the "merger tree" Monte-Carlo algorithm used 
for our calculations, some of the key tests carried out to 
check its validity, and the implementation of two illustra- 
tive models in which only a small fraction of galaxies at 
high redshift harbor a SMBH. The results for the evolution 
of the SMBH occupation fraction and the merger event 
rates as a function of redshift are presented in §3. In §4, 
we discuss some consequences and important limitations 
of our results. 

2. THE MERGER TREE ALGORITHM 
2.1. Dark Matter Halos 

We follow the evolution with redshift of a Press- 
Schechter mass distribution (Press & Schechter 1974; here- 
after PS) of virialized halos at z = 0, using the merger tree 
algorithm described by Sheth & Lemson (1999). For the 
adopted cosmology (described below), we found that this 
algorithm was more accurate in reproducing a PS mass dis- 
tribution after many successive merger steps toward high 
redshift than the algorithm of Somerville & Kolatt (1999). 
In addition, Sheth & Lemson (1999) showed that this algo- 
rithm accurately reproduces the higher order moments of 
the progenitor mass distribution at high redshifts seen in 
numerical simulations. For an assumed PS mass distribu- 
tion of virialized objects at redshift z = 0, the merger tree 
is constructed by picking progenitor masses from an appro- 
priate probability distribution at the next redshift z + Az 
(following the Extended Press-Schechter formalism, here- 
after EPS; Lacey & Cole 1993, 1994), where Az = 0.05 is 
the redshift step of the merger tree. Throughout this work, 
we adopt a ACDM cosmology with Oo = 0.3, fib = 0.04, 
n A = 0.7, ftioo = 0.65, (see, e.g. Bahcall et al. 1999), 
and a mass power spectrum with cr 8 = 0.9 and n = 1 
(see Eisenstein & Hu 1999). We expect our qualitative 
conclusions to remain valid in different cosmologies. Note 
that the h^ scaling for masses and the ft.j~ 3 scaling for 
volumes are not shown explicitly below. 

We have carried out several tests of the merger tree al- 
gorithm, including comparisons to the PS mass function at 
various redshifts. We find that the merger tree algorithm 
reproduces a PS mass function to an accuracy better than 
0.1 dcx till z = 2.5, 0.2 dcx in the range 2.5 < z < 4, and 
0.3 dex in the range 4 < z < 5.0. The discrepancy between 
the merger tree predictions and the EPS predictions for 
the mass function of halo progenitors grows toward higher 
redshifts, and exceeds a factor of 2 at z > 5. Likewise, 
the halo merger rates obtained from the tree agree with 
the EPS rates to within a factor of two for < z < 5. 
A convergence test shows that the redshift step used for 
the calculations (Az = 0.05) does not significantly affect 
the merger rate predictions: the discrepancy between trees 
with Az = 0.01 and Az = 0.05 is largest (~ 25%) at z = 5 
and reduces to nearly zero at z = 0. These tests guarantee 
that the mass functions and merger rates produced by the 
merger tree are accurate to a factor of two or better in 
the redshift range of interest (z < 5). We chose not to ex- 
tend our results beyond z = 5, where the accuracy of the 
merger tree degrades rapidly. Further details and tests of 
the implementation of the merger tree will be given in a 
separate publication (Haiman, Menou & Narayanan 2001, 
in preparation). 

The merger tree we use describes 10 4 halos picked from 



the PS mass function at z = 0, which are subsequently 
broken up into « 8.8 x 10 4 halos by z — 5. Because of 
numerical constraints, the merger tree tracks the merger 
history of halos over a finite range in halos mass, which 
varies with redshift. The maximum mass at any redshift 
corresponds to the most massive halo among the initial 10 4 
halos drawn from a PS distribution at z = (correspond- 
ing to M max w 8.4 x 10 12 M Q ). For the minimum mass, we 
adopt the redshift-dependent value (e.g. Navarro, Frcnk 
& White 1997) 
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corresponding to a virial temperature of the halo of T v ; r = 
10 4 K. This value represents the minimum halo mass below 
which, in a metal-free gas, baryons cannot cool efficiently 
within the age of the universe at that redshift, and there- 
fore cannot form a SMBH. Although this is a reasonable 
choice, the precise value of the limiting cooling mass is un- 
certain (it depends on the mctallicity and H2 chemistry; 
see, e.g. Haiman, Abel & Rees 2000). 

Once a realization of the merger tree has been gener- 
ated, it can be used backwards (i.e., from high z to low 
z) to follow the merger history of dark matter halos which 
potentially contain black holes. We estimate the statisti- 
cal errors on the various quantities calculated below from 
the dispersion in the corresponding quantities among ten 
independent realizations of the merger tree with the same 
parameters. We find that the statistical uncertainties are 
generally small, and are below the expected level of sys- 
tematic errors in the merger-tree method. The merger tree 
predictions can be translated into physical units using the 
equivalent comoving volume that contains the initial 10 4 
halos described by the merger tree. The PS mass function 
predicts 0.579 halos per Mpc 3 at z = 0, in the appro- 
priate mass range M m ; n (z = 0)-M max . Thus, we infer 
that our merger tree simulates a fixed comoving volume of 
AV = 10 4 /0.579 = 1.73 x 10 4 Mpc 3 . 

2.2. Black Hole Seeds And Evolution 

The various models that have been proposed so far for 
the formation of SMBHs in individual galactic nuclei (Rees 
1984 and references therein; Quinlan & Shapiro 1990; 
Eisenstein & Loeb 1995; Silk & Rees 1998; Selwood & 
Moore 1999; Ostriker 2000) generally lack a robust pre- 
diction regarding the presence or absence of a SMBH in 
a galaxy. Here, we explore the predictions for two differ- 
ent illustrative models, in which SMBHs are assumed to 
be present in only a small fraction of all the halos that 
can potentially harbor a SMBH at high redshift [z = 5, 
the initial redshift in the merger tree). We demonstrate 
that: (a) these two widely different models can nonetheless 
be consistent with the near ubiquity of central SMBHs in 
luminous galaxies at z = 0, and (b) the detection of low- 
frequency gravitational waves expected from the merger of 
two SMBHs when the corresponding parent halos merge 
provides a good tool to constrain the extent and proper- 
ties of the population of SMBHs. Throughout this paper, 
we assume that every halo that is modeled by the merger 
tree corresponds to a unique galaxy. 
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Our models are fully specified by the initial "SMBH oc- 
cupation fraction", fumiz, M^), defined as the fraction of 
halos of mass Mh at the initial redshift z (= 5 in our 
merger tree) that harbor a SMBH. In model I, we assume 
that SMBHs are randomly seeded in halos at z = 5 with 
an initial SMBH occupation fraction /em = 3 x 1CP 2 (in- 
dependent of the halo mass above the cooling threshold): 



/BHi(z = 5, M h ) 



3 x 1(T 2 




Mh > M min 
M h < M min 



(z = 5) 
(z = 5), 



( 2 ) 

Here M m j n (z = 5) « 2 x 10 8 M is the cooling mass 
threshold at z = 5. 

In model II, the SMBH occupation fraction, averaged 
over all halos, is /sm = 3 x 10~ 2 as well, but we assume 
that only the most massive halos present at z = 5 are 
seeded with a SMBH 3 : 



/bhiO 



,Mh 



> M crit 
< M cri t, 



(3) 



Here M cr it ~ 5 x 10 9 M is the critical mass such that the 
overall occupation fraction, averaged over all halos with 
masses above M m i n (z = 5), is 3 x 10~ 2 . 

We follow the evolution with cosmic time of this initial 
SMBH population, under the following assumptions: (a) a 
halo merger remnant contains a SMBH only if at least one 
of its progenitors already contained a SMBH prior to the 
merger, and (b) SMBHs in a merger remnant coalesce on a 
timescale that is small compared to the tree redshift step 
(Az = 0.05). We discuss the validity of these assumptions 
in detail in §3. 

In addition, we need to specify the evolution of halos 
whose masses are close to the limiting mass -A^/min ■ As 
mentioned above, we do not track the merger history of 
halos whose mass falls below the limit M m ; n (z) given in 
equation ([I]). As a result, when the merger tree is followed 
from high z to low z, these halos gradually "appear" in 
the tree. In model I, these halos are assigned the same 
seed fraction /em = 3 x 10~ 2 , so that they have the same 
mass-independent probability of containing a SMBH (irre- 
spective of redshift). On the other hand, model II assumes 
that only the most massive halos at z — 5 contain SMBHs. 
Since the newly appearing halos at < z < 5 have small 
masses by construction in model II, we assume that they 
do not contain a SMBH. 

Note that the masses of SMBHs are left unspecified in 
our models. This omission, which greatly simplifies the 
problem, is motivated by the fact that LISA can, in prin- 
ciple, detect the gravitational waves emitted in any merger 
involving black holes with individual masses > lO 3 Af up 
to very high redshifts (see § 3.3.1 below). The accretion 
history of SMBHs, and therefore the evolution of their 
masses, will be explored in a future paper (Haiman, Menou 
& Narayanan 2001, in preparation). In the present work, 
we focus only on the presence or absence of a SMBH in 
any given halo, as a function of redshift. 

3. RESULTS 



3.1. Black Hole Occupation Fraction 

Figures la and lb show the evolution with redshift of 
the black hole occupation fraction /bh in model I and 
model II, respectively. The three panels in each figure 
show, from top to bottom, the black hole occupation frac- 
tion for halos with masses M h > lO n M , M h > 10 10 M Q , 
and for all the halos described by the merger tree (cor- 
responding to the mass range M min (z) < Mh < M max ). 
In both models, the value of /bh, averaged over all halos 
above the limiting mass, is /bhi = 3 x 10~ 2 at z = 5. 
However, the difference in the initial conditions in the two 
models (Eqns. ^ and [|) strongly influences the subsequent 
evolution of /bh via the merging process. Overall, /bh 
computed over the entire mass range M m - m (z) < Mh < 
M max increases with cosmic time because of the numer- 
ous mergers experienced by the halo population. How- 
ever, the evolution of /bh depends strongly on the halo 
mass. For example, in model II, only the most massive 
halos (with M h ^ 5 x 10 9 M Q ) are seeded with a SMBH at 
z = 5, so that /bh for halos with masses Mh > 1O 1O M 
actually decreases with cosmic time because many lower 
mass halos that lack SMBHs become more massive than 
1O 1O M by successive mergers. This decrease in /bh for 
Mh > 1O 1O M is opposite to the evolution seen in model I, 
where the successive mergers of small halos, some of which 
contain SMBHs, increase the value of /bh with cosmic 
time for halos with masses Mh > 1O 1O M . 

Another important effect of the mass-scale dependence 
of the merging process is seen in the value of /bh for the 
most massive halos. In both model I and II, halos with 
masses > 10 11 M Q (corresponding roughly to the masses 
of galaxies brighter than 0.1L* at z = 0) have /bh k5 0.8 
at z = 0. This large SMBH occupation fraction at z = 
can be understood by observing that massive halos possess 
a larger number of early progenitors, and hence the mas- 
sive halos at z = very often had a progenitor seeded with 
a SMBH (in model I). In addition, massive halos are also 
more likely than smaller mass halos to have relatively large 
mass progenitors, and hence are less likely to be "contami- 
nated" by numerous small mass halos without SMBHs (in 
model II). 

We also constructed models with values of the param- 
eter / B Hi = 10 -1 , 10~ 2 and 10~ 3 at z = 5. As might be 
expected, for /BHi = 10 -1 , the resulting values of /bh at 
z = are systematically larger, and nearly unity for the 
most massive (Mh > lO n M ) halos in models I and II. For 
/BHi = 10 -3 , however, the resulting occupation fraction 
/bh at z = for halos of mass M h > 10 12 M Q is only k, 40% 
and w 60% in models I and II, respectively (and even less 
for halos of mass M h > 10 n M o ). For f Bm = 10~ 2 , the 
resulting occupation fraction /bh at z = for halos of 
mass Mh > 1O 11 M is only w 50% in both models. Given 
that the study of Magorrian et al. (1998) probes several 
galaxies down to luminosities w 1O 1O L (roughly corre- 
sponding to a total mass of 1O 11 M for a mass to light 
ratio of 10), 4 the models with /bhi = 10~ 2 do not satisfy 
the observational constraints. Conversely, scenarios with 
a black hole occupation fraction /bh ~ 3 x 10~ 2 or more 



3 We do not consider a model in which SMBHs are preferentially formed in small mass halos because their shallow potential wells do not 
favor such a process (see, e.g., Haehnelt et al. 1998; Silk & Rees 



our models. They show, however, that SMBHs are present in at least a small fraction of lower mass galaxies. 
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at z = 5 are consistent with the observational evidence for 
the presence of a SMBH in nearly all luminous galaxies at 

3 = 0. 

One can ask whether the models are consistent with 
observational constraints from high redshift quasars (see 
also §4 for a discussion of possible additional constraints 
on the models). At z ~ 3, the epoch of maximum activ- 
ity for bright optical quasars, only a fraction ~ 10~ 3 of all 
bright galaxies shows quasar activity (e.g. Richstone et al. 
1998). This extra constraint is not violated in the mod- 
els discussed above. Similarly, the discovery by Barger et 
al. (2001; see also Mushotzky et al. 2000) that a fraction 
~ 10% of all bulge-dominated optically luminous galax- 
ies at z < 3 show hard X-ray activity does not rule out 
the models presented. In fact, we find a large value of 
/bh <; 0.7 for luminous galaxies at z ~ 3 in both classes 
of models considered (see Fig. [j]). This result follows from 
our assumption that SMBHs form mostly (in model I) or 
exclusively (in model II) at high redshift z > 5, and that a 
large fraction £ 80% of halos with M h > 10 11 M Q at z = 
must harbor SMBHs. In order to satisfy the constraint at 
z = 0, at least ~ 3% of halos near the cooling threshold 
at z = 5 must have SMBHs; in model I, the numerous 
mergers between halos in the redshift range 3 < z < 5 
ensures that a fair fraction of all luminous galaxies must 
harbor a SMBH by z ~ 3. This is an important result 
that shows that a fair fraction of all luminous galaxies 
must harbor a SMBH by z ~ 3 if SMBHs preferentially 
form at rcdshifts z > 5. The rapid decline of quasar ac- 
tivity beyond redshifts 2^3 (e.g. Fan et al. 2001) also 
allows for scenarios in which SMBHs preferentially form 
between 3 < z < 5. We did not explore this class of mod- 
els here because the predictions would crucially depend 
on the assumed distribution of formation redshifts for the 
SMBHs. We note, however, that in a scenario with an 
epoch of formation peaking around z <; 3, we find that it 
is possible to have a value of /bh as low as ~ 5 x 10 -2 at 
z = 3 and still satisfy the constraints at 3 = 0. This type 
of scenarios, which would probably lead to reduced merger 
event rates between 3 < z < 5, are better explored with a 
simultaneous attempt at reproducing the observed evolu- 
tion of the quasar luminosity function (Haiman, Menou & 
Narayanan, in preparation). 

The possibility of a relatively rare population of SMBHs 
at high redshifts should be taken into account in studies of 
the evolution of the quasar luminosity function. In terms 
of a fit to a luminosity function, reducing the fraction of 
galaxies that host SMBHs and go through a quasar phase 
can be compensated by either decreasing the mass ratio 
between the SMBH and its host halo by a similar fac- 
tor, and/or assuming a longer activity phase (i.e. a larger 
duty cycle) . Studies of the clustering properties of quasars 
will also have to account for this possibility: the cluster- 
ing would be reduced in scenarios with /bh •C 1 relative 
to scenarios with /bh ~ 1 at high redshift (Martini & 
Weinberg 2000; Haiman & Hui 2000). The problem could 
also be complicated by the fact that in scenarios with rare 
SMBHs at high redshift, the clustering study will have to 
take into account not only the small value of /bh but its 
mass-dependence as well. These issues will be discussed 
further in a future paper (Haiman, Menou & Narayanan 

5 We have checked that the same model without the low— mass halos 
identical merging rates. 



2001, in preparation). 

The general trend shown by the models at z = (due 
to the mass-scale dependence of the merging process) cor- 
responds to larger values of /bh for larger halo masses. 
The transition from halos mostly populated with SMBHs 
to halos without SMBHs occurs over a finite mass range, 
which is characteristic of the extent of the population of 
SMBHs at high redshift. For example, in both model I 
and II, we find that this transition, defined by /bh = 0.5 
(at z = OL occurs at the halo mass Mh ~ 3 x 10 10 M Q (see 
also Fig. [I]), which is approximately three times less than 
the masses extensively probed by the study of Magorrian 
et al. (1998). This property suggests that one way of 
constraining the properties of this population is to search 
for SMBHs in nearby dwarf galaxies with masses smaller 
than have been studied so far (Kormendy & Richstone 
1995; Magorrian et al. 1998). 

3.2. Rates of Black Hole Mergers 

In this section, we show that the detection of gravita- 
tional waves emitted when two SMBHs merge is an ef- 
ficient technique for determining the extent of the pop- 
ulation of SMBHs over a wide range of masses and red- 
shifts. It is possible to calculate the cosmic rate of halo 
mergers containing SMBHs from the merger tree. If the 
corresponding SMBHs coalesce (see §3.3.2), these merger 
events are potentially detectable by a low-frequency gravi- 
tational wave experiment such as LISA. The total number 
of events A A occurring in the tree during one redshift step 
(A3 = 0.05) corresponds to AN / Az / AV events per Mpc 3 
per unit redshift (used here as a time coordinate). Given 
the adopted cosmology, this event rate can be translated 
into the number of detectable merger events per unit time 
per unit redshift (used here as a space coordinate; see e.g., 
Kolb & Turner 1990; Haehnelt 1994 for similar deriva- 
tions): 

dN AN , \ 2 j2 / \ 
dzdt-A^AV^ 1 + Z ^ d ^ < 4 > 

where cLa{z) is the angular diameter distance and c is the 
speed of light. 

Figure 0b shows the event rates in model I and model II. 
We also show in Figure ||a the event rates expected in a 
third "maximal" model in which we assume that every 
single halo described by the merger tree harbors a SMBH 
(i.e., /sHi = /bh = 1 at any redshift, independent of halo 
mass). The predicted event rates differ both in terms of 
the total number of events and in their distribution with 
redshift. 

The maximal model shows that the merging rate of all 
the halos described by the merger tree decreases steeply 
with cosmic time (Fig. ||a). For the adopted cooling mass 
(Eq. 0), the merger tree predicts as much as ~ 550 merger 
events per year per unit redshift at z — 5 and less than 
~ 30 merger events per year per unit redshift below 3 = 1. 
In model I, on the other hand, mergers of halos containing 
SMBHs are much rarer because of the reduced population 
of SMBHs in halos (/bhi = 3 x 10 -2 , independent of mass). 
However, since mergers gradually increase the value of /bh 
with cosmic time (Fig. nk), the probability of a merger in- 
volving two SMBHs initially increases with redshift. 5 The 

'appearing" from below our mass resolution limit predicts essentially 
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black hole merger event rate peaks at z ~ 2.5 where S5 3 
events per year per unit rcdshift are predicted (Fig. |2|b; 
dotted line). This is the most pessimistic of the models 
considered here, with a total event rate ~ 9 yr _1 out to 
z = 5. 

In model II, by construction, mergers of halos with 
SMBHs initially occur only for halos more massive than 
w 5 x 1O 9 M . As halos merge together, this limit is in- 
creased further at z < 5, so that mergers occur on a typ- 
ical mass scale which is growing with cosmic time. Halo 
merger rates generically peak at a characteristic redshift 
that varies inversely with the halo mass scale (i.e., the 
lower the halo mass, the higher the peak redshift); it peaks 
at a redshift beyond 5 for ~ 1O 1O M0 halos (Lacey & Cole 
1993; see also Figure 5 in Haiman & Menou 2000). The 
combination of these two properties leads to a decrease 
in the merger rate with cosmic time in model II (Fig. |^b; 
dashed line), where the typical halo mass involved in a 
merger with SMBHs increases with cosmic time. The 
merger tree predicts ~ 9 black hole merger events per year 
per unit redshift at z = 5, decreasing to less than ~ 0.5 
merger events per year per unit redshift below z = 1. A 
total of ~ 15 events per year out to z — 5 is predicted. 

3.2.1. Low- Frequency Gravitational Waves 

SMBHs are prime targets for LISA (Folkncr 1998) as 
sources of low-frequency (1 — 10~ 4 Hz) gravitational waves 
when they coalesce (see Press & Thorne 1972; Thorne 1987 
and references therein) . The merger tree allows us to make 
more definite estimates of the black hole coalescence rates 
than has been previously done (see Thorne & Braginsky 
1976; Haehnelt 1994; 1998). These predictions must take 
into account the spectral coverage and sensitivity of the 
instrument. However, following Haehnelt (1994; see also 
Folkner 1998), we note that LISA should be able to de- 
tect the gravitational signal from two coalescing SMBHs 
up to very high redshifts (z > 10) provided they both have 
masses <; 1O 3 M . 

The masses of SMBHs are left unspecified in the merger 
tree implementation described in §2. A simple ansatz for 
the SMBH masses is to assume that they follow, up to 
z — 5, the relation determined locally: 



M BH = 1-2 x 10 S M Q 



3.75 



(5) 



.200 km s- 1 - 

where o~ e is the velocity dispersion of the stellar bulge 
(Gebhardt et al. 2000). The relation between the stellar 
and dark matter velocity dispersions can be specified based 



on the observed radial stellar distributions (p s ta 



oc r 



and assuming an isothermal DM halo (pdm oc r ), which 
yields a e = (3/2) -1 / 2 ctdm (see, e.g., Turner, Ostriker 
& Gott 1984). Adopting this relation, we find that all 
the SMBHs described by the merger tree have masses 
> 10 3 Mq (the minimum SMBH mass described by the 
tree at z = 5 is « 1.3 x 1O 3 M , according to Eqns. |Q 
and H). In this case, all the merger events counted in 
Figure g are detectable by LISA. 

If instead of equation (0) , we adopt the steeper relation 
of Merritt & Ferrarese (2001a) 



M BH = 1-3 x 10 8 M e 



( 5» ) 

V200 km s" 1 / 



4.72 



(6) 



we find that our smallest halo harbors a SMBH with a 
mass <J 100 M at Z = 5, below LISA's detection limit. 
In this case, the rate of events detectable by LISA could 
be substantially reduced. We note, however, that both 
equations (JsJ) and (|^) have been obtained in nearby galax- 
ies. In extending these relations to higher redshifts, one 
needs to ensure that the total amount of mass in SMBHs 
is consistent with the observed luminosity of the quasar 
population (Soltan 1982). This implies that in scenarios 
with rare SMBHs at high redshift, the typical SMBH in 
a halo with a given mass will have to be more massive 
than in scenarios where SMBHs are more widespread. As 
a result, the above relations between velocity dispersion 
and black hole mass derived locally should scale differ- 
ently at high redshift, making SMBHs more massive (and 
more detectable) in our models where /bh < 1- This sug- 
gests that only scenarios in which SMBHs are common at 
high rcdshift will be affected by the above mass/sensitivity 
threshold for the rate of merger events detectable by LISA. 
This issue, which requires a detailed accretion model for 
the SMBHs besides their merger history, will be addressed 
in more detail elsewhere (Haiman, Menou & Narayanan 
2001, in preparation). 

The gravitational wave signal provides a wealth of in- 
formation on the emission source (e.g. Thorne 1987). For 
coalescing black hole binaries, the observables are the bi- 
nary orbital period P OT b (twice the quadrupolar wave pe- 
riod), its time derivative P rb, and the gravitational wave 
amplitude h. Modulo the unknown binary inclination 
and source-detector orientation for individual sources, one 
can put meaningful constraints on the binary mass and 
its distance to the observer via a combination of these 
three observables (e.g. Schutz 1986). Thus, the quantity 
{h P rb / Porb) provides a direct estimate of the distance to 
the source, independent of its mass; the value of h then 
constrains the mass of the system once the distance is 
known. 6 This raises the exciting possibility that LISA 
could not just detect the events predicted in Figure |[ but 
could also discriminate between coalescing SMBHs at high 
and low redshift, constrain the masses involved in the var- 
ious merger events and, all together, probe the character- 
istics of the population of SMBHs as a function of redshift. 



3.2.2. Black Hole Mergers 

The SMBH merger event rates shown in Figure | were 
derived assuming that SMBHs coalesce efficiently, i.e. 
shortly after the merger of their parent halos (§2.2). Our 
current knowledge of the merger process is limited, and 
various plausible physical effects could substantially mod- 
ify its efficiency. 

In their seminal work, Begelman, Blandford & Rees 
(1980) showed, for a specific massive example, that the 
coalescence timescale for two SMBHs in a galactic merger 
product could be very long, i.e., comparable to or larger 
than a Hubble time. Numerical simulations show that 
dynamical friction between the SMBH and the surround- 
ing star cluster, while being an important effect before 



.200 km s- 

6 In a cosmological context, P or b and P OT b are affected by the cosmological redshift and the distance becomes a luminosity distance, but 
equally good constraints can be derived assuming the cosmology is known (e.g. Schutz 1989). 
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the phase of binary evolution dominated by gravitational 
wave emission, is rather inefficient, so that the coalescence 
of the two SMBHs is not an obvious outcome (Rajagopal 
& Romani 1995; Makino 1997; Quinlan & Hernquist 1997; 
Merritt, Cruz & Milosavljevic 2000). However, this conclu- 
sion is complicated both by the possible acceleration of the 
coalescence process due to the "wandering" of the SMBH 
binary in the central stellar cluster — a process that is 
difficult to follow numerically (Quinlan & Hernquist 1997; 
Merritt et al. 2000) — and by the possible influence of an- 
other merger on the preexisting stellar distribution (Roos 
1988; Polnarev & Rees 1994; see also Begelman et al. 1980 
for a discussion of gas infall effects). 

If SMBHs are brought together by successive halo merg- 
ers at a rate higher than the rate at which they can coa- 
lesce, the lowest mass SMBH is likely to be ejected out of 
the nucleus of the merger remnant by the slingshot mecha- 
nism (Saslaw, Valtonen & Aarseth 1974; Hut & Rees 1992; 
Xu & Ostriker 1994). This black hole ejection is more 
probable in the maximal model presented above (where 
/BHi = 1) than in models I and II, where SMBHs are com- 
paratively rare. For example, in the maximal model, a 
halo with mass > 10 12 M Q at z = typically has experi- 
enced about 10 2 to 10 3 mergers between < z < 5 with 
other halos containing a SMBH. In fact, we find that in the 
maximal model, ~ 10% of all the mergers occurring within 
a redshift step Az — 0.05 (corresponding to« 1.6 x 10 7 yrs 
at z — 5 and « 7 x 10 8 yrs at z = for the adopted cos- 
mology) are multiple mergers involving three or four halos 
all of which contain a SMBH (indicating therefore very 
short times between two successive mergers). Note that 
most of the mergers happen at high redshift in this model, 
as shown by Fig. ^a. In models I and II, these "multiple" 
mergers within a redshift step Az = 0.05 are exceedingly 
rare and the cumulative number of mergers of halos with 
SMBHs involved in the building-up of the halos at z = 
is much reduced. For example, halos of mass > 10 11 M Q 
at z = have experienced up to 50 mergers in model I, 
and up to 80 mergers in model II, between < z < 5, 
with other halos containing a SMBH (although the ~ 20% 
of halos that still do not contain a SMBH at z — have 
not experienced any such merger; see Fig. [l]b). In all three 
models, there is therefore the possibility for slingshot ejec- 
tions via three-body interactions, unless the SMBH merg- 
ing process is very efficient (i.e. coalescence time <C than 
the typical time between successive mergers, at any red- 
shift). Analytical arguments and numerical simulations, in 
the context of dark matter halos being made of SMBHs, 
suggest that many of the ejected SMBHs would remain un- 
detectable and bound to their parent galaxy, while some of 
them could escape to infinity (see, e.g., Hut & Rees 1992, 
Xu & Ostriker 1994). This raises the possibility of the 
existence of a population of free-floating SMBHs located 
outside galaxies. 

4. DISCUSSION 

Our description of the merger history of SMBHs is lim- 
ited to redshifts z < 5: we find that the merger tree fails to 
accurately reproduce the PS mass function at higher red- 
shifts (§2.1). In the model with a SMBH in every halo sat- 
isfying our cooling criterion (the maximal model), and in 
model II, the halo merger event rates are strongly peaked 
toward high z, suggesting even more events at z > 5. How- 



ever, this simple extrapolation toward higher redshifts of 
the expectations for LISA is limited by the tendency for 
SMBHs to be less massive at high z (thereby moving the 
coalescence signal out of the frequency range of LISA), 
and by the fact that bound objects do not form beyond 
z ~ 20 — 30 in popular CDM cosmologies (see Barkana & 
Loeb 2001 for a recent review). 

Our models also suffer from the uncertainty on the min- 
imum cooling mass M mm (Eq. jl]]). Because of the steep 
slope of the halo mass function, most of the halos described 
by the merger tree have masses close to M m j n , so that 
the predictions for the merger event rates in the maxi- 
mal model (/bh = 1) and, to a much lesser extent, in 
model I are affected by the value of M m j n (more merger 
events are expected for smaller M m i n ). Nevertheless, we 
expect our qualitative conclusion to hold with different 
values of M m i n : (1) a few percent of the massive halos at 
z = 5 should harbor a SMBH to satisfy the constraints 
at z = and (2) if SMBHs do not preferentially form be- 
tween 3 < z < 5, then a fair fraction of luminous galaxies 
must contain a SMBH at z ~ 3, the epoch of peak quasar 
activity. 

We also note that the models I and II considered above 
are highly idealized. In model II, it is postulated that 
SMBHs form exclusively in high-mass galaxies at z > 5 
and, implicitly, that no new SMBHs form at < z < 5. 
The assumption in model I, namely that a small univer- 
sal fraction of the halos described contain SMBHs, could 
represent the scenario envisioned by Eisenstein & Loeb 
(1995), in which SMBHs form only in halos with unusu- 
ally low specific angular momentum. We note that, in this 
scenario, halos in which SMBHs form should experience 
fewer mergers because of their comparatively poor envi- 
ronments and /eHi should somewhat depend on the halo 
mass. These two effects are not included in our model I. 

As we have shown, the merging rates scale strongly with 
halo mass, so that in a more realistic model the fraction 
/bh would be an increasing function of halo mass at almost 
any redshift (and in particular at the arbitrary starting 
redshift z = 5 chosen here) because of previous mergers 
experienced by the halos at even larger redshifts. It is also 
possible that SMBHs continuously form with cosmic time, 
with some complex dependence on redshift. Nonetheless, 
these idealized models represent two plausible extremes 
of the black hole occupation functions, serving to illus- 
trate the essentials of the merging process and showing 
that SMBHs could be rare in halos at high redshift with- 
out violating constraints at z = 0. Indeed, since the same 
models with an additional, significant rate of SMBH for- 
mation with cosmic time or an initial parameter /em that 
is an increasing function of halo mass would predict larger 
values of /em at z = for the large mass halos, the two 
models explored here are useful extreme cases. 

A source of uncertainty for the models presented here is 
that real galaxies (e.g. those in clusters) can correspond to 
sub-halos of a larger collapsed halo (e.g. the cluster). This 
invalidates the one-to-one correspondence between galax- 
ies and Press-Schechter halos, at least for galaxies resid- 
ing in groups or clusters (see White, Hernquist & Springel 
2001 for recent results on the multiplicity function of sub- 
halos). We note that the sample of galaxies of Magorrian 
et al. (1998), that is used here as the main observational 
constraint for the models, is non-uniform, containing some 
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galaxies from the Virgo cluster, the Local Group, as well as 
other galaxies. Implicit in our work is the assumption that 
nearly all galaxies in the field (those with a one-to-one cor- 
respondence to Press-Schechter halos) harbor SMBHs. If 
future observations reveal that SMBHs are less common in 
field galaxies, this could allow even more extreme models 
than considered here (i.e. with even rarer SMBHs). 

The results of Magorrian et al. (1998) have been quan- 
titatively challenged, in that the SMBH masses are now 
thought to be a factor of a few smaller than originally 
estimated by Magorrian et al. (see, e.g., Merritt & Fer- 
rarese 2001b; Gebhardt et al. 2000). We note that these 
uncertainties do not effect our conclusions here, since in 
our models, we only prescribe the presence or absence of a 
SMBH. Although this has not been established to date, if 
SMBHs were found in the future to be rarer (not just less 
massive) than estimated by Magorrian et al., the case for 
rare SMBHs at high redshift could be made even stronger 
(i.e. the fraction of SMBHs populating galaxies at high 
redshift could be even smaller than considered here). 

Evolution models with rare SMBHs can be further con- 
strained by requiring that they fit the observed luminosity 
function of quasars. Since only a small fraction, /bh "C 1, 
of high-redshift halos are populated with SMBHs in these 
models, those halos which do contain SMBHs had to har- 
bor more massive SMBHs in these models than they would 
have in standard evolution models with /bh = 1. In stan- 
dard models, the mass that constitutes a central SMBH at 
z — is distributed among all the halo progenitors at some 
earlier redshift z, each of whom harbors a SMBH. In com- 
parison, in models with rare SMBHs, the same material 
resides at redshift z in a fraction ~ /bh of the progeni- 
tors, so that the typical ratios of black hole to halo mass 
are larger in models with rare SMBHs by a factor ~ /^. 
In standard models (see, e.g., Haiman & Loeb 1998; Kauff- 
mann & Haehnelt 2000), the SMBH has a mass which is 
a small fraction, ~ 10~ 2 — 10~ 3 , of the total gaseous mass 
of the parent galaxy. This places an approximate limit of 
/bh <; 10~ 2 — 10~ 3 on the maximum rarity of SMBHs. 

Successful evolutionary models are also required to re- 
produce the tight correlation between the SMBH mass and 
the stellar velocity dispersion of the host galaxy bulge ob- 
served locally (Ferrarese & Merritt 2000; Gebhardt et al. 
2000). We emphasize that in the models presented here, 



the masses of high-z SMBHs or their accretion history were 
not specified and that there is considerable freedom in ad- 
justing both of these quantities in an evolutionary model. 
For example, in the models of Kauffmann & Haehnelt 
(2000) for the evolution of active nuclei, the fraction of 
the gas that is accreted onto the SMBH during a merger 
is chosen to scale with redshift. We expect that a suitable 
prescription for the accretion histories and/or the initial 
SMBH masses can lead to a tight A/bh — cr relation as 
observed locally, and will pursue this question in a future 
work. 

The above points demonstrate the need for a better 
characterization of the population of SMBHs as a func- 
tion of redshift. Wc highlighted the possibilities and lim- 
itations of the detection of gravitational waves associated 
with SMBH mergers for this purpose. It would be very use- 
ful to test the extent of the population of SMBH via other 
means. The frequency of SMBH binaries at the centers of 
nearby galaxies may well be different for the various mod- 
els discussed above, although this frequency depends on 
the efficiency with which SMBHs coalesce, and is difficult 
to estimate theoretically. A more promising test could be 
to search for SMBHs at the centers of local dwarf galaxies 
with masses smaller than currently probed, since our mod- 
els predict that nuclear SMBHs may not be ubiquitous in 
galaxies below a certain characteristic mass scale at z = 0. 
In the future, searches for nuclear activity at high redshift 
will be useful in constraining models with rare SMBHs, be- 
cause they provide lower limits to the fraction of galaxies 
harboring SMBHs (some of them being possibly inactive). 
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Fig. 1. — The black hole occupation fraction /bh is shown as a function of redshift in model I (a) and model II (b). In each case, the three 
panels (top to bottom) correspond to halos with mass > 1O 11 M0, > 1O 1O M0 and all the halos (meaning roughly > 1O 8 M0) described by 
the merger tree. The large mass halos (M > 1O 11 M0) nearly all harbor a SMBH at z = despite the low overall value of /bh- Standard 
deviations are indicated by error-bars. 
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Fig. 2. — Predictions for the black hole merger event rates (under the assumption of very efficient black hole coalescence) in scenarios where 
(a) all halos susceptible to harbor a black hole and (b) only a fraction /bh = 10 — 2 of those, do contain a black hole. In (b), the dotted line 
corresponds to model I and the dashed line to model II; the event rates in (b) were smoothed over 8z = 0.2 for better rendering. Standard 
deviations are indicated in each case by the upper and lower thin lines. Various effects which could modify these rates arc discussed in the 
text. 



